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ABSTRACT 

As the quality and quantity of astrophysical data continue to improve, the precision with which certain 
astrophysical events can be timed becomes limited not by the data themselves, but by the manner, standard, and 
uniformity with which time itself is referenced. While some areas of astronomy (most notably pulsar studies) 
have required absolute time stamps with precisions of considerably better than 1 minute for many decades, 
recently new areas have crossed into this regime. In particular, in the exoplanet community, we have found that 
the (typically unspecified) time standards adopted by various groups can differ by as much as a minute. Left 
uncorrected, this ambiguity may be mistaken for transit timing variations and bias eccentricity measurements. 
We argue that, since the commonly-used Julian Date, as well as its heliocentric and barycentric counterparts, 
can be specified in several time standards, it is imperative that their time standards always be reported when 
accuracies of 1 minute are required. We summarize the rationale behind our recommendation to quote the 
site arrival time, in addition to using BJDtdb- the Barycentric Julian Date in the Barycentric Dynamical Time 
standard for any astrophysical event. The BJDtdbIs the most practical absolute time stamp for extra-terrestrial 
phenomena, and is ultimately limited by the properties of the target system. We compile a general summary of 
factors that must be considered in order to achieve timing precisions ranging from 15 minutes to 1 /is. Finally, 
we provide software tools that, in principal, allow one to calculate BJDxdb to a precision of 1 /is for any target 
from anywhere on Earth or from any spacecraft. 
Subject headings: Barycentric, Heliocentric, Julian Date, UTC, TDB, TT, TDT, TAI, light time 



1. INTRODUCTION 

Precise timing of astrophysical events is one of the funda- 
mental tools of astronomy, and is an important component of 
essentially every area of study. There are two basic sources 
of uncertainty in timing: the astrophysical data characterizing 
the event, and the time stamp with which the event is refer- 
enced. Unfortunately, since the accuracy of the time stamp is 
something that is often taken for granted, the improvements 
in data are sometimes not accompanied (or are not uniformly 
accompanied) by the requisite improvements in accuracy of 
the time stamp used. This situation can lead to confusion, or 
even spurious inferences. 

Timing plays a particularly important role in the study of 
exoplanets. Indeed, many of the ways in which exoplanets are 
discovered involve the detection of transient or time-variable 
phenomena, including the radial velocity, transit, microlens- 
ing, and astrometry techniques. Furthermore, in some cases 
much can be learned about planetary systems from the pre- 
cise timing of these phenomena. As examples, the measure- 
ment of terrestrial parallax in microlensing allows one to in- 
fer the mass of the primary le ns and so the planetary com- 
panion (e.g.. iGould et al.l 12009 ), and one can constrain the 
eccentricity of transiting planets by comparing the ti mes of 
primary transits and secondary eclipses (e.g ., Demin g et alj 
120051; ICharbonneau et alll2005i; iKnutson et al...2007.) . Possi- 
bly the most promising application of timing in exoplanets, 
however, comes from transit timing variations (TTVs). With 
an exquisitely periodic phenomena like transiting planets, we 
will be able to measure many effects using the departures from 
strict periodicity, s uch as the gravitationa l perturbations from 
additi o nal planets (lMiralda-Escudell2002t iHolman & Murray! 
loos'; 'Agol et all l2005h . troians ([Ford & Gaudi' 120061 
Ford &Holman" l2007h. and moons (IK ipping 2009), stel- 
lar quadrupoles jMiralda-Escudd I2002D . tidal deforma- 



tions, genera l relativistic precession (iJordan & BakosI 120081: 
iPaL&Kocsis' 2008; Hevl & Gladm anI |2()07D. orb ital decay 



(e.g..[ Sasselov .2003). and proper motiori dRafikovll2009l) . 

Because of their great potential, TTVs have become the fo- 
cus of many groups. The typical data-limited transit timing 
precisions of most observations are around 1 minute, with 
the best transit time precision yet achieved of a few seconds 
(Pont et al. 2007). However, as discussed above, accuracies of 
transit times are limited not only by the data themselves, but 
also by the time stamp used. In order to make these difficult 
measurements useful, it is critical that a time stamp be used 
that is considerably more accurate than the uncertainty due to 
the data themselves. Furthermore, since thorough characteri- 
zation of TTVs will require the use of all available data span- 
ning many years from several groups, this time stamp must 
be stable in the long term, and all groups must clearly convey 
how it was calculated. 

Unfortunately, we have discovered that, in the exoplanet 
community, the Julian Date (JD) and its geocentric (GJD), 
heliocentric (HJD), and barycentric (BJD) counterparts are 
currently being quoted in several different, and often unspec- 
ified, time standards. In addition, the site arrival time and 
its time standard is not quoted. This general lack of homo- 
geneity and specificity leaves quoted time stamps ambiguous 
at the 1 minute level. More alarmingly, the most commonly- 
used time standard, the Coordinated Universal Time (UTC), is 
discontinuous and drifts with the addition of each leap second 
roughly each year. 

The pulsar community has solved the problem of precise 
timing well beyond the level that is currently necessary for ex- 
oplanet studies, and we can benefit from the techniques they 
have developed over the past 40 years. In particular, their 
current state of the art program (TEMP0 2) models pulsar ar- 
rival t imes to 1 ns precision (Hobbs et al.l2006l;lEdwards et alj 
120061) . This program is highly specialized and generally can- 
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not be applied outside of pulsar timing observations, but many 
of the effects they consider are relevant to optical observers in 
the exoplanet community. 

In this article, we summarize the effects one must consider 
in order to achieve timing accuracy of 1 /is - well beyond the 
accuracy that will likely be required by the exoplanet commu- 
nity for the foreseeable future. 

Section |2] provides the background required to understand 
each of the effects that could change the arrival time of a pho- 
ton. They are listed in order of decreasing magnitude, so latter 
subsections can be ignored for low-precision measurements. 

Section |3] discusses the practical limitations to achieving 
high-precision timing. We begin with the effects which may 
cause errors that are comparable to or exceed the BJD cor- 
rection. These should be read and understood by everyone. 
We continue with remaining effects, in order of decreasing 
magnitude, which can be ignored for low -precision (> 30 ms) 
measurements. We conclude § 3 by listing additional effects, 
the errors due to which are negligible (< 1/is). 

We begin § |4] by detailing the procedure one must follow 
in order to calculate the BJDtdb, which is designed to be a 
useful reference for those already familiar with the concepts 
of precision timing. In the latter part of this section, we de- 
scribe our particular IDL and web-based implementation of 
this procedure. 

Lastly, in the Appendix, we discuss some of our specific 
findings about the time stamps currently in use and how these 
are calculated throughout the exoplanet community. 

While we focus on the effects of timing on the opti- 
cal/infrared exoplanet community, timing precision of order 
1 minute is necessary for many oth er areas, such as th e study 
of rapidly rotating white dwarfs tEuchner et al.ll2006l) . This 
article should be equally applicable in such cases. 

2. THEORETICAL TEVIING PRECISION 

The biggest source of confusion comes from the fact that 
time standards and reference frames are independent from one 
another, even though there are many overlapping concepts be- 
tween the two. We will use the following terminology: "ref- 
erence frame" will refer to the geometric location from which 
one could measure time - different reference frames differ by 
the light-travel time between them; "time standard" will refer 
to the way a particular clock ticks and its arbitrary zero point, 
as defined by international standards; and "time stamp" is the 
combination of the two, and determines the timing accuracy 
of the event. 

The BJDjDB, the time stamp we advocate, can be calculated 
using the equation: 
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Fig. 1 . — Difference between BJDtdb and the uncoiTected JDtt (see text 
for definitions) over the course of a year. We plot the correction for a 
maximally-affected object on the ecliptic for an observer at latitude of zero 
degrees. We exclude all points where the object has an airmass greater than 
three and the Sun is higher than -12° in order to highlight observing biases. 



done accurately enough for their purpose, the site arrival time 
(e.g., the JDuTc) should always be quoted in addition to the 

BJDxDB- 

2.1. Reference Frames: The R0mer Delay 

Due to the finite speed of light, as the Earth travels in its 
orbit, light from an astrophysical object may arrive early or be 
delayed by as much as 8.3 minutes from the intrinsic time of 
the extraterrestrial event. This is called the R0mer delay, Ar, 
in honor of Ole R0mer's demonstration that the speed of light 
is finite. Since most observers cannot observe during daylight, 
a bias is introduced and in practice the delay (as distinct from 
the early arrival time) is only as much as 7 minutes, for a peak- 
to-peak variation of 15 minutes. Figure [T] shows an example 
of this effect for a maximally affected object on the ecliptic. In 
order to show the observational bias, our example assumes the 
object is at O'' right ascension and 0° declination. This curve 
shifts in phase with ecliptic longitude and in amplitude with 
ecliptic latitude. We also place our observer at the Earth's 
equator, but note that the asymmetry will be larger at different 
latitudes. 

The solution to this problem is to calculate the time when 
a photon would have arrived at an inertial reference frame. 
This time delay is the dot product of the unit vector from the 
observer to the object, n, and the vector from the origin of the 
new reference frame to the observer, r 



BJD 



TDB '■ 



JDuTc + Arp, + Ar + Asp, + A, 



■S0- 



^£0: 



(1) 



where JDutc is the Julian Date in Coordin ated Universal 
Time; Ar0 is the R0mer Dela y, dis cussed in § 12.11 Ac is the 
clock correcti on di scussed in § 12.21 Asq is the Shapiro delay 
discussed in § 12.31 and Aeq is the Einstein delay, discussed 
in §1231 

The order of these terms is such that they are of decreas- 
ing magnitude, so one need only keep the terms up to the 
precision required. The timing precision required by current 
exoplanet studies (^ 1 s) requires only the terms up to and 
including Ac. 

Because future Solar System ephemerides may enable more 
precise calculations of the arrival time at the Barycenter, or in 
order to allow others to check that the original conversion was 



Arq = 



(2) 



where c is the speed of light and n can be written in terms of 
its right ascension (a) and declination (S), 



(cos((5)cos(a)^ 
cos((5) sin(a) 
sm{6) J 



(3) 



This equation is general as long as a, 8, and r are in the 
same coordinate system (e.g.. Earth mean equator J2000) and 
the object located at (a, S) is infinitely far away. Other 
forms of this equation in the literature assume that we have 
the angular coordinates of the new origin or that the Earth 
and the new origin are in the same plane (e.g., .Binnendij fl 
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Fig. 2. — Difference between the BJDtdb and the HJDtdb for a 
maximally-affected object on the ecliptic. The primary periodicity is due 
to Jupiter and the secondary periodicity is due to Saturn. 



[T960HHend"en & Kaitchu ckll98alHirs"hfeld & Sin nott'iggT), 
but we explain in §|4]why this form is most practical for cal- 
culating the delay. 

The HJD, which uses the Sun as the origin of the new refer- 
ence frame, is only accurate to 8 s because of the acceleration 
of the Sun due primarily to Jupiter and Saturn (Fig. |2]). It 
was popular when people first began considering this effect 
because it is relatively simple to calculate from tables without 
a computer (Landolt & Blondeau 1972) . and remains popu- 
lar because self-contained algori thms exist to approx imate it 
without any external tables (e.g.. lDuffett-Smithll 19891) . How- 
ever, because the HJD is not useful when accuracies of better 
than 8 s are needed, most of the algorithms in use today use 
approximations that are only precise at the 1 s level, and it 
becomes impossible to back out the original JD from the HJD 
unless we know the exact algorithm used. 

Because of these problems, the HJD was formally depre- 
cated by International Astronomical Union (lAU) Resolution 
A4 in 1991, in favor of the BJD, a time referenced to the Solar 
System Barycenter (SSB). 

The analogous correction to the R0mer delay in our Solar 
System can also be significant in the target system. We refer to 
this as A«. For example, for transiting planets with a « 0.06 
AU, Ar can be as large as 30 s. In general, the position of 
the planet during primary transit has become the unspoken 
standard reference frame for transiting planets, while the host 
star's photosphere is the unspoken standard for Radial Veloc- 
ity (RV) planets. 

In theory, the timing would be much more stable in the tar- 
get's barycentric reference frame, but the accuracy with which 
we can convert to this frame depends on the measurements of 
the system. Since different observers may use different values 
as measurements improve, quoting the JD in the frame of the 
target's barycenter may obfuscate the long term reliability of 
timing. Therefore, we argue it is better to quote Julian Date in 
the SSB reference frame, and correct for Ar only when com- 
paring observations at different phases in the planet's orbit. 

This Ar correction is not necessary for TTVs of the primary 
transit, since the planet is always in the same phase. Neverthe- 
less, we should explicitly state the object's reference frame to 
avoid any potential ambiguity, particularly when comparing 
any combination of primary transits, secondary transits, RVs, 
and another primary transit of a different planet in the same 



system, when it may not be obvious which origin is being 
used. 

For RV measurements, which are taken at many different 
phases, the effect is much smaller and can generally be ig- 
nored because the star's orbit around the barycenter is small. 
For a typical hot Jupiter, (i.e., a Jupiter mass planet in a 3 day 
orbit around a solar mass star), the maximum time difference 
in the RV signal (for an edge-on orbit) is 20 ms, which would 
change the measured RV by ^ 50 fim s~' . While planets far- 
ther out will cause a larger timing offset, the difference in the 
measured RV is even smaller 

2.2. Time Standards: The clock correction 

To be clear, the J D can be specified in many time standards 
(ISeidelmannll 1 992h . and while the lAU has made no explicit 
statement regarding the allowed time standards of the GJD, 
HJD, or BJD, their meaning in any given time standard is un- 
ambiguous. Unfortunately, they have been specified in many 
standards, usually implicitly. 

However, the particular time standard used affects how use- 
ful the time stamp is as an absolute reference. We must be 
careful not to directly compare BJDs or HJDs in different 
time standards, as each has different offsets, periodic terms, 
and/or rates, which can introduce systematic errors of over 1 
minute. For this reason, it is critical that any stated BJD or 
HJD also specify the time standard used when one-minute ac- 
curacies are important, and the uncertainty of a time that is 
quoted without a standard should be assumed to be at least 1 
minute. 

First, it may be useful to summarize the relevant standards 
of time: 

Universal Time, UTl: - Defined by the mean solar day, and 
so drifts forward and backward with the speeding and 
slowing of the Earth's rotation. Generally, it slows due 
to the tidal braking of the Moon, though changes in the 
Earth's moment of inertia and complex tidal interac- 
tions make its exact behavior unpredictable. It is rarely 
used directly in astronomy as a time reference, but we 
mention it for context. 

International Atomic Time, TAI: - Based on an average of 
atomic clocks all corrected to run at the rate at the geoid 
at K, with 1 s equal to "9,192,631,770 periods of 
the radiation corresponding to the transition between 
the two hyperfine levels of the ground state of the cae- 
sium 133 atom," as defined by Resolution 1 of the thir- 
teenth meeting of the Conference Generale des Poids et 
Mesures (CGPM) in 1967. This definition is based on 
the duration of the Ephemeris Time second, which was 
previously defined as 1/3 1,556,925.9747 of the tropical 
year for 1900 January at 12 hours Ephemeris Time 
by Resolution 9 of the eleventh CGPM in 1960. TAI 
is the fundamental basis for many other time standards, 
and is the default time standard of the Sloan Digital Sky 
Survey. 

Coordinated Universal Time, UTC: - Runs at the same 
rate as TAI, except that it is not allowed to differ from 
UTl by more than 0.9 s. Every 6 months, at the end of 
3 1 December and 30 June, the International Earth Rota- 
tion and Reference Systems Service (lERSflmay elect 

' At |http://www.iers.org| 
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to add (or subtract) a leap second to UTC in order to 
keep it within 0.9 s of UTl. UTC is therefore discontin- 
uous and drifts relative to TAI with the addition of each 
leap second, which occur roughly once per year. As of 
January 2009, the current number of leap seconds, A^, 
is 34. The full table of leap seconds is available online 
and is typically updated several months in advance of 
when an additional leap second is to be addecQ. UTC is 
the current international standard for broadcasting time. 
As a result, when a modern, network-connected com- 
puter's clock is synchronized to a Network Time Pro- 
tocol (NTP) server, it will be in UTC. Thus, this is the 
system of time most familiar to astronomers and non- 
astronomers alike (modulo time zones and daylight sav- 
ings time). 

Universal Time, UT: - An imprecise term, and could mean 
UTl, UTC, or any of several other variations. In gen- 
eral, such imprecise language should be avoided, as the 
potential ambiguity is up to 1 s. In the context of a time 
stamp, it is likely UTC, but some people may intention- 
ally use UT to imply 1 s accuracy. While explicitness 
is preferred (i.e., UTC ± 1 s), any time stamp quoted in 
UT should be assumed to be uncertain at the 1 s level 
unless the time standard has been independently veri- 
fied. 

Terrestrial Time, TT(TAI): - A simple offset from TAI of 
32.184 s released in real time from atomic clocks and 
never altered. This offset is to maintain continuity be- 
tween it and its predecessor, the Ephemeris Time (ET). 

Terrestrial Time, TT(BIPM): - A more precise version of 
TT(TAI). The International Bureau of Weights and 
Measures (BIPM) reanalyzes TT(TAI) and computes 
a more precise scale to be used for the most demand- 
ing timing applications. The current difference between 
TT(TAI) and TT(BIPM) is - 30/is, and must be inter- 
polated from a table maintained by the BIPM and pub- 



lytl 



lished online with a 1 month dela 



Terrestrial Time, TT: - Sometimes called Terrestrial Dy- 
namical Time (TDT), can refer to either TT(TAI) or 
TT(BIPM). From this point on, we will not make the 
distinction, but when accuracies of better than 30 /is 
are required, TT(BIPM) must be used. 

Barycentric Dynamical Time, TDB: - Corrects TT for the 
Einstein delay to the geocenter, A^q, which is the de- 
lay due to time dilation and gravitational redshift from 
the motions of the Sun and other bodies in the Solar 
System. The conversion from TT to TDB cannot be 
written analytically, but is usually expressed as a high- 
order series approximation (Irwin & Fukushima 1999). 
The difference is a predominantly a periodic correction 
with a peak-to-peak amplitude of 3.4 ms and a period 
of I yr. TDB was slightly modified by lAU Resolu- 
tion B3 in 2006, converging on the same definition as 
the JPL Ephemeris Time, Tgph, also called Coordinate 
Time (CT) in the JPL ephemerides of Solar System ob- 
jects. 

^ At'ftp://maia.usno.navy.miiyser7/tai-utc.dat' 
3 At ftp://tai.bipm.org/TFG/TT %)28BIPM%29. 



Barycentric Coordinate Time, TCB: - Physically and 
mathe matically equ ivalent to the TDB as defined in 
2006 ( IStandishlll998i ). and differs only by an offset and 
rate of about 0.5 s yr~' due primarily to time dilation in 
the Sun's gravitational potential. TDB and TCB were 
roughly equal at 1977 Jan 1.0 TAI, and now differ by 
about 16 s. 

Caution must be always be exercised, however, as these def- 
initions are subject to change at any point, though usually with 
a few year's notice. 

Assuming the time is measured according to current defi- 
nition of UTC, then the clock correction, from UTC to TDB, 
can be written as the sum of the corrections from UTC to TAI, 
TAI to TTQ, and TT to TDB: 



Ac = N+32AS4s + iTDB-TT). 



(4) 



Of course, if one wishes to express the BJD in another 
time standard (or start with something other than UTC), the 
clock correction would change accordingly. However, not ev- 
ery time standard is well-suited to precise, astrophysical time 
stamps, and the use of any time standard other than TDB 
should be viewed simply as an adequate approximation to 
TDB. 

Most readily available programs that calculate the time 
stamp assume that the user has already applied the UTC-to- 
TT part of this clock correction, which is often not true. We 
feel this assumption has contributed the widespread confusion 
regarding time stamps. 

As this last point is our primary motivation for writing this 
article, we elaborate here on the effects of time standards on 
the reliability of time stamps. For the sake of simplicity, we 
only discuss the effects of time standards on the BJD. Each of 
these effects also applies to the HJD, though the improvement 
in accuracy of the time stamp is negligible compared to the 
accuracy of the HJD reference frame for all but the UTC time 
standard. 

The least preferred, though most commonly used, time 
standard for the BJD is UTC (BJDutc), and is equivalent to 
ignoring Ac altogether. Because UTC is discontinuous and 
drifts with the addition of each leap second, comparing two 
BJDuTC time stamps could result in spurious differences if 
any leap seconds have been introduced between observations. 
Therefore, 1 s timing accuracies cannot be achieved using the 
BJDuTC over a span that straddles the addition of one or more 
leap seconds (roughly 1 yr). Figure [3] shows the difference 
between the BJDutc time stamp with the uniform BJDtdb 
time stamp (described below) from 1961 Jan 1, when UTC 
was defined (though its definition has evolved over the years), 
to 2010 December 31, the furthest future date for which the 
value of UTC can be accurately predicted at this writing. 

BJD in TT (BJDtt), which is equivalent to ignoring the 
(TDB-TT) term in equation (|4|, corrects for the discontinuity 
and drift introduced by leap seconds and is appropriate for 
timing accurate to 3.4 ms. 

BJD in TDB (BJDtdb) is usually the best time stamp to use 
in practice, as it further corrects the BJDtt for all known ef- 
fects on the motions, and therefore rates, of our atomic clocks. 
While BJDtdb is not perfect, any more accurate time stamp is 
unique for each target. 

BJD in TCB (BJDtcb) corrects for the gravitational po- 
tential, primarily from the Sun, which causes the clock to 

^ Ignoring the TT(BIPM) correction. 
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Fig. 3. — Difference between the uniform BJDtdb and the BJDutc- It 
shows the discontinuities and slow drift in BJDutc due to the addition of leap 
seconds. Without connecting for these, relative timing between two reported 
values of BJDutc can only be trusted over short time scales. 

run slower than it otherwise would. However, if one is con- 
cerned about effects of this magnitude, the analogous cor- 
rection of putting it in the gravitational potential of each ob- 
served object is also required. Since these two rates are small 
(10~**), and opposite in sign, we believe it is best to ignore 
BJDtcb except perhaps as an intermediate step in calculat- 
ing the target-specific frame. Technically, the use of TCB was 
recommended by lAU Resolution B 1.3 in 2000. However, be- 
cause of the greater practicality of using TDB (see § |4]i, and 
the drifting difference between TCB and TDB and TT, we 
believe its use will only lead to confusion without any fore- 
seeable benefit to the exoplanet community. 

2.3. Shapiro Delay 

The Shapiro delay, A^ ( IShapirolll964l) . is a general rela- 
tivistic effect in which light passing near a massive object is 
delayed. For an object at an angle 6 from the center of the 
Sun, the Shapiro delay is 



A 



so ■ 



2GMq 

5 — log( 1 - cos fc*) . 



(5) 



This can be as large as 0. 1 ms for observations at the limb 
of the Sun, but for objects more than 30° from the Sun, the 
correction is less than 20 /is. 

There is also the analogous correction. As, for the target 
system. Similar to Ar, A^ depends on the measurements of 
the target, which may be refined over time. Therefore, the 
time should generally be quoted without As, but include it 
when comparing times where this could be significant. 

2.4. Einstein Delay 



As we discussed in § 12.21 relativity dictates that the motion 
of the observer influences the rate at which the observed clock 
ticks. The use of TDB corrects for an observer moving with 
the geocenter, but in reality we observe from the surface of 
the Earth or from a satellite, for which there is an additional 
term to Aeq: 



A£0 = 



ro-V(s 



(6) 



Fig. 4. — Difference between the BJDtdb calculated for an object at R.A. = 
O*" and an object at R.A. = O*" 0™ 4^* observed at the same time. This difference 
is as large as 200 s for a 10° offset. 



there is an analogous correction, Ag, for the target system, 
which should be ignored when quoting the time but included 
when comparing times if necessary. 

3. PRACTICAL CONSIDERATIONS 

Of course, the accuracy of the output time stamp is only 
as good as our assumptions and the inputs needed to specify 
the time standard and reference frame. Here we discuss their 
effects on the accuracy of the time stamps. The first four sub- 
sections, through § 13.41 have no reasonable upper bound. In 
each case, the accuracy of the inputs must be evaluated de- 
pending on the accuracy of time stamp required. The later 
three subsections through § 13.71 are organized in decreasing 
magnitude, all of which ca n be ignored for accuracies no bet- 
ter than 21.3 ms. Finally, § l3.8l discusses effects of less than 1 
/is that we have ignored. 

3.1. Coordinates 

An error of 1 ' in the position of the target amounts to a tim- 
ing error of as much as 0.28 s in the BJD (Fig. |4|i. Such error 
would be common if the coordinates of the field center are 
used instead of the specific object's coordinates. In particular, 
if doing a survey, one may wish to assign the same BJD to all 
objects in a given frame. However, with a 10° offset, which is 
possible with some wide-field transit surveys, the error can be 
as large as 200 s for objects at the edge of the field. An error 
of 0.25" will yield 1 ms timing offsets, and 0.25 mas accuracy 
is necessary for 1 /is timing. 

3.2. Computer Clock 

The accuracy of a typical computer clock depends on its in- 
trinsic stability, the computer workload, its operating system, 
and the reliability of the network connection. Older comput- 
ers with a parallel port CCD interface may produce unreliable 
timing because the clock may slow or stop completely during 
CCD readout. Without any special effort, a modern Windows 
machine with a network connection is accurate to '^ 2 sec- 
ondfl and with third-party software like Dimension 4E we 
have found it to be stable to 0.1 s. An NTP-synchronized 
Linux machine is typically accurate to ^ 50 ms. 



Here, r^ is the location of the observer with respect to the 
geocenter, and Vq is the velocity of the geocenter. Again, 



At http://support.microsoft.com/kb/939322j 
At http://www.thinkman.com/dimension4j 
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Of course, the stability of the clock only sets a lower limit 
on the absolute accuracy of the time recorded in the FITS im- 
age header. NTP synchronization attempts to measure and 
compensate for network latency, but the accuracy of time 
stamps also depends on the particular software package tak- 
ing the image and the hardware it uses, which is difficult to 
calibrate. Unless independently verified, the time recorded 
in image headers should not be trusted to better than 0.25 s. 
However, various solutions exist for higher precision timing, 
such as GPS -triggered shutters. 

In particular, it is worth emphasizing the 1 s error in the 
Hubble Space Telescope (HST) clock and potential 6.5 s er- 
ror in the Kepler clock, described in more detail in the Ap- 
pendix, both of which have already achie ved data-limited 
transit timing precisions of that order (e.g.. iPont et alj|2"007t : 
iKippi ng & Bakos 2010). 

3.3. Flux-weighted mean time of exposure 

When calculating the time of exposure of an image, we typ- 
ically use the time at midexposure. However, the precise time 
of exposure that should be used is the flux-weighted mean 
time of exposure. The magnitude of this error depends on the 
intrinsic stability of source, the stability of the atmosphere, 
and the exposure time. In the diabolical case of a cloud com- 
pletely covering the object during one half of the exposure, 
the error could be as large as half the exposure time. For 
a typical Hot Jupiter that dims by 1% over the course of a 
15 minute ingress, the error introduced into the time stamp 
during ingress or egress by using the mid exposure time is 
0.25% of the exposure time - 150 ms for a 1 minute expo- 
sure. Near the peak o f some high-magnification microlensing 
events (i.e.. iGouldet al. 2009), the flux may double in as little 
as 6 minutes. During such instances, using the mid exposure 
time will result in a time stamp error of as much as 1/4 of the 
exposure time -15 s for a 1 minute exposure. 

3.4. Plane-Wave Approximation 

Equation ^ assumes the object is infinitely far away, and 
therefore the incoming wavefronts are plane waves. In real- 
ity, the wavefronts are spherical, which introduces a distance- 
dependent, systematic error. The maximum error introduced 
by the plane-wave approximation is 1000 s for the MoorQ, 
100 s for the Main Asteroid Belt, 5 s for the Kuiper Belt, 1 ms 
at the distance of Proxima Centauri, and 150 ns at the distance 
to the Galactic Center 

The fully precise equation, assuming spherical wavefronts, 
is 



\r'+nd\ —d 
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where d is the distance from the observer to the target. In 
such instances, the distances must be derived from precise 
ephemerides for both the target and observer 

Although this form is generally applicable, it is not gen- 
erally practical because at large distances, double precision 
floating point arithmetic cannot reliably recover the small dif- 
ference between |r + nii| (the distance from the bary center to 
the target) and d. 

' The maximum error for the Moon is not because of the departure from 
the plane wave, but because the plane-wave formalism can place the Moon 
on the wrong side of the SSB. 
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Hour Angle 

Fig. 5. — Difference between the BJDtdb calculated at the geocenter and 
at the precise location of the observer on the surface of the Earth. While 
geometrically, this effect will oscillate between ±21.3 ms with a period of 1 
sidereal day, we exclude points when the Sun is above -12° and object is at 
z > 3, which introduces a large observational bias. 



To solve this problem, the pulsar community (e.g.. lLorimeri 
20080 uses the two-term Taylor expansion of equation (|7]l 
about \/d = Q: 
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One may recognize the first term as plane wave approximation 
(eq. ||2l)- In practice, the accuracy of this two-term approxi- 
mation exceeds the accuracy of the "exact" calculation using 
double precision at a distance of 10,000 AU (0.05 pc). At the 
distance of Proxima Centauri, the accuracy of this approxima- 
tion is at worst 1 ns. 

However, this Taylor expansion is divergent when li < |rj. 
It should never be used for objects inside 1 AU and may still 
be inadequate for other objects inside the Solar System. It has 
a maximum error of 1 day for the Moon, 20 s for the Main 
Asteroid Belt, and 40 ms for the Kuiper Belt. In these cases, 
using the exact formula (eq. Q) may be easiest. 

Therefore, we recommend that for precise calculations of 
any Solar System body, equation ^ should be used. For bet- 
ter than 1 ms timing of any object outside the Solar System, 
equation (O should be used. 

3.5. Geocenter 

Most readily available time stamp calculators use the posi- 
tion of the geocenter, rather than the location of the observer 
on the surface of the Earth. Neglecting the light-travel time 
from the surface of the Earth to the center introduces a 2 1 .3 ms 
amplitude variation with a period of 1 sidereal day. In prac- 
tice, most observers can only observe their targets at night, 
creating a systematic bias of between 8 ms and 21.3 ms (Fig. 

3.6. BJD'jjjfj versus BJDjjtc 
Usually, when people calculate the BJDutc, they neglect 
Ac, and input JDutc for algorithms designed to take JDtt- 
This effectively uses the positions of the Earth, Sun, and plan- 
ets offset in time by 32.184 -i- N s to calculate the correction. 

** Note: the sign of the second term in Lorimer 2008, eq. 9, should be - 
rather than +. 
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Fig. 6. — Difference between the correctly calculated BJDutc an'l the com- 
monly calculated BJDyjj, using the positions of the Earth, planets, and Sun 
delayed by 32.184 + N s. It shows that the difference can safely be ignored 
for ~ 15 ms precision, and therefore the approximate BJDtt or BJDtdb can 
be recovered from cun'ently published BJDutc simply by adding 32. 1 84 + N 



When the BJDutc is calculated in this manner, we denote it as 
BJDjjTc- The correct way to calculate the BJDutc would be 
to first calculate the BJDtdb, then subtract the Ac correctioifl 
Figure |6] shows, for a maximally affected object on the eclip- 
tic, an example of the difference between the BJDjj.j,p and the 
fully correct BJDutc- Fortunately, this amounts to at most a 
13.4 ms difference (though growing with the UTC-TT differ- 
ence), which is below the precision of most clocks and the 
geocentric correction that is usually ignored. Therefore, to an 
accuracy of ~ 50 ms, one can safely say BJDtdb ~ BJDtt ~ 
BJDuTc+Ac, making it easy to convert currently published 
values of BJDutc '^^ '^^e superior BJDtdb- 

3.7. Computer Precision 

Representing JD as a double precision floating-point num- 
ber limits the accuracy to about 1 ms, and any operation done 
on the full JDs will be even less accurate. Many programs 
require the use of a Reduced or Modified Julian Date, and/or 
can return the JD to BJD offset in seconds, but care must be 
taken at every step of the way never to store the full JD as a 
double precision number if 1 ms precision is required. 

3.8. Negligible (< l^is) effects 

The Shapiro delay occurs for other bodies as well, but ob- 
servations at the limb of Jupiter only delay light by 200 ns. 

Typical modern, commercial GPS units use the World 
Geodetic System (WGS84), which is referenced to the Inter- 
national Terrestrial Reference System (ITRS) with an error of 
about 15 m, which amounts to a 50 ns error in the time stamp. 

The index of refraction of the atmosphere is not exactly 1 
and changes with its composition, temperature, and pressure, 
which changes the speed of light. However, the largest rea- 
sonable deviation due to this effect is only tens of ns. 

The pulsar community must specify a frequency-dependent 
dispersion measure delay. At radio wavelengths (21 cm), the 
delay can be as much as 1 s, but the dispersion delay con- 
tributes at less than 1 ^s shortward of ^ 300 /im. 

' However, the BJDutc is a poor approximation to BJDtdb and should 
not be used. 



4. CALCULATING THE BJDtdb 

The most practical way to precisely calculate the BJDtdb 
time stamp is using JPL's DE405 ephemerio- It contains 
the position of thousands of bodies in the Solar System, in- 
cluding the Sun, planets, spacecraft, moons, asteroids, and 
comets. It is oriented to the International Celestial Reference 
Frame (ICRF), which is consistent with the FK5 system at 
J2000.0 within the 50 mas error of FK5, and has its origin 
at the S SB with its axes fixed with respect to extragalactic 
objects (I Arias et al. 1995). Therefore, it is recommended to 
use the 3-space Cartesian coordinates retrieved from the JPL 
DE405 ephemeris directly with the J2000 object coordinates 
in equation (|2]l. 

The following is an outline of the steps required to properly 
calculate the BJDtdb - 

1. Calculate the midexposure time in JDutc- Most FITS 
image headers give DATE-OB S in UTC at the begin- 
ning of exposure. If high precision is required (< 1 s), 
read the caveats about clock precision in § l3.2l carefullv: 
depending on the sky conditions or variability of the 
object, one may need to account for the flux-weighted 
mean time of exposure. 

2. Convert the midexposure time to JDtdb by applying Ac 
(eq. ID). For times accurate to 3.4 ms, one can use the 
simpler JDtt and calculate BJDtt (using the positions 
of planets delayed by the TT-TDB offset). The differ- 
ence between BJDtt ^i^d BJDtt is no more than ^ 200 
ns, which is well below the precision of the BJDtt- If 
better than 30 /is precision is required, the TT(BIPM) - 
TT(TAI) offset must be applied. 

3. Retrieve the JPL ephemeris of the observing station for 
the times spanning the observing window. JPL's HORI- 
ZONS svsterrPl is designed for this. To return inputs 
for use with equation (|2) and J2000 target coordinates: 
select "Vector Table", the SSB as the coordinate ori- 
gin, and "Earth mean equator and equinox of reference 
epoch" for the reference plane. This will return the 
Cartesian coordinates of the observing station with re- 
spect to the SSB in the J2000 Earth mean equator ref- 
erence frame at steps as small as 1 minute in CT, which 
is the same as TDB. 

4. If the observing station is on Earth and better than 20 
ms timing is required, another ephemeris must be gen- 
erated (from HORIZONS) for the observer's position 
with respect to the non-rotating geocenter, and added to 
the geocentric positions. We note that the precise con- 
version from latitude, longitude, and elevation to the 
Cartesian coordinates with respect to the non-rotating 
geocenter is not trivial, and requires tables of measured 
precession and nutation of the Earth. 

5. Interpolate the positions of the observing station to each 
midexposure time JDtdb - 

6- Input the interpolated X, Y, and Z positions of the ob- 
serving station, and the target's J2000 Earth mean equa- 
tor coordinates into equation (|2|i. Depending on the dis- 
tance to the target and the precision required, equation 

'" At' ftp://ssd.jpl.nasa.gov/pub/eph/e xport/DE405/de405iom.ps I 
" At http://ssd.jpl.nasa.gov/7horizons] 
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(|7]) (and the target's ephemeris) or equation (O may be 
required. One must be careful to use sufficiently accu- 
rate target coordinates. 

7. If greater than 0.1 ms precision is required, apply the 
Shapiro correction (equation ||5l). An ephemeris of the 
Sun is required, which can be generated from HORI- 
ZONS. 

8. If greater than 1 /is precision is required, apply the ad- 
ditional Einstein correction for the observing station's 
position with respect to the geocenter The geocentric 
velocity is required, which can also be given by HORI- 
ZONS. 

4.1. Our Code 

Our IDL code implementing this procedure is available on- 
lin43 It requires the JDutc at midexposure and target coor- 
dinates (a, S) in J2000 as inputs. We outline its procedure 
here. More explicit details, as well as the calling procedure 
and dependencies, are commented inside the code. 

We compute Ac using Craig Markwardt's TAI_UTC pro- 
gran{3 to read the leap second table, and his TDB2TDT 
program to compute the TT-TDB correction, which uses a 
791 -term Fairhead and Bretagnon analytical approximation 
to the full numerical integr ation, with an error of 23 ns 
(iFairhead & Bretagnonll 19901) . 

Our code will automatically update its leap second table the 
first time it runs after every January 1 or July 1, but this re- 
quires a periodic Internet connection and the use of the wget 
program. It will terminate on failure to update, but this protec- 
tion can be bypassed for those that elect to (or have to) update 
their table by hand. 

By default, we ignore the - 30fis TT(BIPM) - TT(TAI) cor- 
rection, which would require a constant Internet connection, 
would not apply to data acquired in the previous month, and 
is likely negligible for most applications. However, our code 
can optionally correct for it if an up-to-date file is supplied. 

To read and interpolate the ephemeris from JPL, 
we use Craig Markwardt's routines JPLEPHREAD and 
JPLEPHINTERP for the Earth, Sun, and other planets. If the 
observing station is space-borne, the smaller ephemeris used 
with those programs does not include satellites, so we use an 
Expect script to automate a telnet session to the HORIZONS 
system and automatically retrieve the ephemeris, which we 
quadratically interpolate to the desired times using IDE's 
INTERPOL. The accuracy of this interpolation depends on 
how quickly the object's position is changing and the step size 
of the ephemeris. HORIZONS can only return '- 90,000 data 
points per query, so the smallest step size (1 minute) limits 
the calculation to a range of 60 days. For the geocenter, a 100 
minute step size is sufficient for 60 ns accuracy, but, for exam- 
ple, a 2 minute step size is required for 1 /is accuracy for the 
HST (though it is still limited by its clock accuracy). We have 
found that a 10 minute step size is adequate for 1 ms timing 
for most objects and allows a range of nearly 2 yr. 

If the observer is on the Earth, and the coordinates (lati- 
tude, longitude, and elevation) are given, we correct for the 
additional delay. If no observer-specific information is given, 
we assume the observer is at the geocenter, and the result will 
be biased by ~ 10 ms (Fig. |5]l. 

'^ At http ://astroutils . astronomy. ohio-state . edu/time ' 

'^ At http://cow.physics.wisc.edu/ 7o7Ecraigm/idl/ephem.htmI. 



If the target's ephemeris can be returned by HORIZONS 
and its unique name is given, we use our Expect script to gen- 
erate its ephemeris too, and calculate the exact Ar0 (eq. ||7l). 
If not, and instead the distance is given, we use the two-term 
approximation to the spherical wave solution (eq. IS)). Oth- 
erwise, we use the plane wave approximation (eq. J!)). 

Lastly, we include the Shapiro correction and the additional 
Einstein correction due to the position of the observer with 
respect to the geocenter, either from the surface of the Earth 
(if given the coordinates), or the spacecraft. 

In the geocentric case, our code agrees with BARYCENnto 
200 ns (peak to peak) and the authors of BARYCEN report 
that their code agrees with FXBAR^iQto 1 /is. The ephemeris 
we generate for a location on the surface of the Earth agrees 
with HORIZONS to 20 nano-lt-s, and the geocentric BJDs we 
calculate from HORIZONS ephemeris agree with the BJDs 
we calculate using Craig Markwardt's routines within 10 ns. 

The near-exact agreements between these methods is not 
surprising, and do not necessarily indicate that they are accu- 
rate to better than 1 /is. Our code was inspired by BARYCEN 
and both rely on Craig Markwardt's routines (the difference 
comes from the fact that we index the JPL ephemeris with 
JDtdb instead of JDtt), and all methods use JPL's DE405 
ephemeris. 

The primary advantage of our code is that it includes the 
JDuTC to JDtt correction (but can optionally ignore it). The 
choice of starting with JDutc is a departure from what is typ- 
ically done with such time stamp calculators, but we feel this 
is a far more robust starting point. The current confusion has 
shown that many assume JDuxc as the starting point, which 
is likely due to a lack of explicitness in the programs and/or 
unfamiliarity with various time standards. Our hope is that 
people are unlikely to make the opposite mistake (assume the 
input should be JDtt instead of JDutc) since our code is very 
explicit and calculating the JDutc is almost always a trivial 
calculation from the DATE-OBS FITS header keyword. 

Additionally, our code can easily correct for the observer's 
position on the Earth or from a spacecraft, and can include the 
spherical wave correction. 

In order to schedule observations, even ±10 minute pre- 
cision is generally good enough, and one can approximate 
B JDtdb ~ JDutc; for more demanding observing schedules, 
we provide software to iteratively calculate the reverse cor- 
rection. 

Along with the IDL source code, we provide a web-based 
interface to our codeo. though not every feature is enabled. 
Specifically, it is limited to 1 ms precision, can only do one 
target at a time, only does the plane wave approximation, and 
is limited to 10,000 JDs at a time. Those with applications 
for which these features are too limited should download our 
source code and run it locally. 

5. CONCLUSION 

Timing of transient events is a powerful tool for characteriz- 
ing many astronomical phenomena. In the field of exoplanets 
in particular, the search for variations in the times of primary 
transits and secondary eclipses, or transit timing variations, is 
one of the most promising new techniques for studying plan- 
etary systems. 

The accuracy with which transit times, and indeed any tran- 



At http://astro.uni-tuebingen.de/software/idl/aitlib/astro/barycen.pro] 



'^ At http://heasarc.nasa. gov/lheasoft/ftools/fhelp/fxbary.txt 

'* At http://astroutils.astronomy.ohio-state.edu/time/utc2bjd.htnill 
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sient phenomenon, can be measured is limited not only by 
the data themselves, but by the time stamp to which the tran- 
sit time is referenced. As the quality of transit timing data 
crosses the threshold of 1 minute precision, the precise time 
standard and reference frame in which event times are quoted 
becomes important. Achieving uniform and accurate time 
stamps with accuracies of better than 1 minute that can re- 
liably be compared to one another requires extraordinary care 
in both our techniques and our terminology. We have found 
that the time standards adopted by various groups that mea- 
sure transit times can differ by as much as a minute, and are 
typically left unspecified. As these ambiguities can be sig- 
nificant compared to the timing precisions that are quoted, 
they may therefore lead to spurious detections of transit tim- 
ing variations or biased eccentricity measurements. 

Here we have summarized the effects one must consider in 
order to achieve timing precision of 1 /is. We argue that the 
BJDtdb is nearly the ideal time stamp, being as reliable as 
any time stamp can be without being unique to each target 
system. On the other hand, BJDutc and the HJD in any form 
should be avoided whenever possible. Most importantly, we 
emphasize that the time standard should always be explicitly 
stated. Any time stamp that is quoted without a time standard 
should be assumed to be uncertain to at least 1 minute. Unless 



the time standards used in programs or algorithms have been 
independently confirmed, one should avoid using ones that do 
not precisely specify the input and output time standard. 

In addition, the arrival time at the observing site along with 
its time standard (e.g., JDuxc) should also be specified. This 
will remove any ambiguity in the time stamp, allow others to 
apply improved corrections should more precise ephemerides 
become available in the future, and allow others to check that 
original conversion was done accurately enough for their pur- 
pose. 

Finally, we have written an IDL program for general use 
that facilitates the use of BJDtdb to an accuracy of 1 /^s, pro- 
vided that the inputs are sufficiently precise, and we provide 
a web-based interface to its most useful features. 
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APPENDIX 

We looked in detail at several readily available tools for the BJD/HJD calculation, and have been in contact with many people in 
the exoplanet community and the help desks for several major space telescopes. We summarize our findings here to demonstrate 
how easily errors of up to 1 minute can be introduced and to stress the importance of specifying the time precisely. We caution 
the reader not to trust our general findings for specific cases, but always to confirm what has been done in each case where 1 
minute timing accuracy is required but the time standard has not been specified explicitly. 

Software 

Ground-based observers have typically used one of the following methods to calculate the BJD. However, most FITS image 
headers give the DATE-OBS and TIME-OBS keywords in UTC. We have found that most people, when starting with JDutc end 
up quoting HJD[jjf, or BJDyj^,. 

HORIZONS 

JPL's HORIZONS ephemeris calculator, which is used by many to calculate the BJDs from space telescopes, and can be used to 
calculate ground-based BJDs, returns the time in JDct= Teph= JDjoswhen the ephemeris type is "Vector Table". Any conversion 
that uses a HORIZONS ephemeris in JDtdb but indexes it with JDutc, as had been done by several people we spoke with, will 
calculate BJD(j.j.(-, which can be offset from the true BJDutc by up to 13.4 ms (as shown in Fig. |6]l, and offset from the uniform 
BJDtdb by more than 1 minute (as shown in Fig. |3]l. 

IRAF 

IRAF's set j d calculates the HJD, but calls for UT, which is likely to be interpreted as JDutc- In this case, it would calculate 
the drifting quantity HJDjj^j.^. If TT were used instead, it would calculate the HJDtt, accurate to ^8 s. 

IDL 

The IDL routines HELIO_JD (for HJDQ from the IDL Astronomy Library curated by Wayne Landsman, and BARYCEN (for 
BJD), from the Institut fiir Astronomic und Astrophysik IDL Library, maintained by Jorn Wilms, both call for the GJD, which, we 
remind the reader, can be specified in any time standard. Often, this is interpreted as JDutc, in which case they would calculate 
HJD(jjf, or BJDfjjQ, respectively. If TT were used, they would calculate the HJDtt, accurate to ~8 s, or BJDtdb, accurate to the 
geocentric correction (21.3 ms). 

AXBARY 

NASA's High Energy Astrophysics Science Archive Research Center (HEASARC) created the tools FXBARY and later the 
improved version FAXBARY, both of which call AXBARY to calculate the BJD. Their documentation is precise and correct, but 
quite long and may be difficult for the uninitiated to follow. Therefore, it would not be surprising for users of AXBARY to input 
either UTC or TT, in which case they could generate either the BJDfj-j.^ to the accuracy of the leap seconds or the BJDtdb to the 
accuracy of the geocentric correction (21.3 ms). 



At |http://idlas tro.gsfc.nasa.gov/ftp/pro/astro/helioJd.prQ I 
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Online Tools 

Currently, common Google results turn up various applets, spreadsheets, programs, or algorithms to calculate the HJD that 
explicitly call for JDutc or JDut as an input. Unless explicitly mentioned otherwise, it is usually safe to assume the time 
standard used as input will be the time standard used throughout their calculation. Thus, these algorithms and applets will very 
likely calculate HJDjjjf;,. However, they are perfectly capable of calculating HJDtt if given JDtt as an input. 

Space telescopes 
EPOXI 

EPOXI has the mid-exposure time BJDtdb in the header for the intended pointing under the FITS header keyword KPKSSBJT. 
This can be used directly, as long as one is careful about the intended target, so it would be very surprising if a BJD from EPOXI 
was not BJDtdb- 

We recalculated the BJDtdb using the HORIZONS ephemeris as described in ^ and an example FITS header given to us by 
the EPOXI help desk. In their example FITS header, they pointed at the Moon, which is not infinitely far away, so we must use 
equation |7] for Arq. With this method, we agree with the KPKSSBJT header value to '^ 1 ms, the limit of the precision of the 
keyword. 

We also redid the BJDtdb calculation of HAT-P-4b as described in the Report on the Calibration of EPOXI spacecraft timing 
and reduction to Barycentric Julian Date of August 2009 by Hewagama et al.. We find agreement with the quoted KPKSSBJT 
FITS header key word to 47 ms. While this is much better than the 0.41 s difference calculated by Hewagama et al. (a difference 
they attribute to "cumulative rounding limits" in their method), we believe our method to be far more precise. However, we 
could not obtain access to the original headers and were unable to determine the source of the discrepancy. Given the very good 
agreement with the calculation of the Moon above, our best guess is that the target coordinates used by the EPOXI pipeline 
differed from the published values for HAT-P-4b. The 47 ms difference could be explained by a ^ 27" discrepancy in R.A., a 
^ 45" discrepancy in declination, or some combination thereof. 

Chandra X-Ray Observatory 

Chandra stores their DATE-OBS keyword in TT. Their more precise TSTART and TSTOP keywords are expressed in seconds 
after 1998 Jan 1, 00:00:00 TT. This departs from what is typically done, which may lead to confusion, but it makes the conversion 
to a uniform time stamp much more straightforward and less likely to drift by the leap seconds. They provide extensive directions 
online to calculate the BJDtdb using AXBARY, so it is likely that anyone using Chandra who quotes a BJD is using BJDtdb- 

Hubble Space Telescope 

The FITS headers of HST state that their DATE-OBS and TIME-OBS keywords are UT. We contacted the i?5rhelp desk for 
clarification, since UT is ambiguous. The //5rhelp desk response stated that their clock reports UTC accurate to ~10 ms, but 
"due to variabilities and quantization in the particular science instruments' operations, the actual time light begins falling on the 
detector is not known to better than about ^1 second, ± 50% (rough estimate)." It is possible for the HST engineering team 
to calibrate these variations, but they have limited resources and have no plans to do so. It is thought that this error is some 
combination of random and systematic errors, but the precise breakdown is unknown. 

This potential ^1 s systematic error may have important implications for the reliability of the transit times quoted with HST 
observations, most importantly, the 3 s error of the transit time of HD 189733b (Pont et al. 2007). 

HST does not calculate the HJD or BJD at any point, leaving the calculation up to each individual observer. Our experience 
with ground based observers suggests that most people will end up quoting an HJDy^f, or BJDyj^;;. 

Kepler 

The Kepler Data Release Notes 2 describe how to calculate the BJD from UTC, but do not include the correction to TDB. They 
mention the HORIZONS ephemeris, but neglect to mention that its output time is in CT, not UTC; thus it appears they calculate 
BJD(j-j.(^, though we were unable to confirm this. In addition, the time stamp uncertainty may be much larger than typical, so it is 
worth quoting from the Kepler Data Release Notes 5 (released 2010 June 4): The advice of the DAWG [Data Analysis Working 
Group] is not to consider as scientifically significant relative timing variations less than the read time (0.5 s) or absolute timing 
accuracy better than one frame time (6.5 s) until such time as the stability and accuracy of time stamps can be documented to 
near the theoretical limit. 

Spitzer Space Telescope 

The Spitzer pipeline calculates the HJDutc for the intended pointing (presumably the target) at the end of the exposure, 
subtracts the full exposure time, and records the result in the header as HJD. Depending on the exposure time, this will produce 
roughly a 10 ms effect similar to that shown in Figure|6l and depending on how close the intended pointing was to the object of 
interest, may produce a ~ 0.1 s effect similar to FigureH) However, this effect is negligible compared to both the ^8 s accuracy 
of the HJD (Fig. |2]) and the number of leap seconds that may have elapsed between observations (Fig. O. 

One typically quotes the HJD at the midexposure time. Since Spitzer quotes the HJDutc at the beginning of the exposure, 
using the unmodified Spitzer HJDs would produce a systematic offset of half the exposure time, though experienced observers 
correct for this. 

Also contributing to this confusion, the FITS header keyword UTCS_OBS is incorrectly documented. While the documentation 
states that it is seconds after J2000 ET, it is actually seconds after January 1st, 2000 12:00 UTC H- A^ - 32. Therefore, trusting the 
documentation as is will unwittingly lead to a difference of A^+32.184 s. 
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Most of the people we have asked opted to calculate their own BJD using the HORIZONS ephemeris. However, they have 
typically quoted the BJD(j.j.(;.. 
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